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The diffusion quantum Monte Carlo method is extended to solve the old the- 
oretical physics problem of many-electron atoms and ions in intense magnetic 
fields. The feature of our approach is the use of adiabatic approximation wave 
functions augmented by a Jastrow factor as guiding functions to initialize the 
quantum Monte Carlo prodecure. We calculate the ground state energies of 
atoms and ions with nuclear charges from Z = 2, 3, 4, . . . , 26 for magnetic field 
strengths relevant for neutron stars. 

Keywords: Quantum Monte Carlo; Atoms; Neutron Stars; Magnetic Fields. 



1. Introduction 

The discovery of features in the X-ray spectra of the thermal emission spec- 
tra of the isolated neutron star IE 1207^'^ and three other isolated neutron 
stars has revived the interest in studies of medium-Z elements in strong 
magnetic fields. The reason is that the observed features could be due to 
atomic transitions in elements that are fusion products of the progenitor 
star. However, to calculate synthetic spectra for model atmospheres, and 
thus to be in a position to draw reliable conclusions from observed spec- 
tra to the elemental composition of the atmosphere and the distribution 
of elements on different ionization stages, accurate atomic data for these 
elements at very strong magnetic fields 10^ to 10^ T) are indispensible. 

While the atomic properties of hydrogen and, partly, helium at such field 
strengths have been clarified in the literature over the last 25 years (for a 
detailed list of references see, e. g., Ref. 3), for elements with nuclear charges 
Z > 2 only fragmentary atomic data exist with an accuracy necessary for 
the calculations of synthetic spectra. 
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Wc have tackled"^ the problem by adapting the diffusion Monte-Carlo 
method (DQMC)^"^ to the case of neutron star magnetic fields. This 
method has the advantage that ground-state energies can be determined 
practically free from approximations. 

2. DQMC for neutron star magnetic fields 

The basic idea of DQMC is to identify the ground state wave function 
$0 (-R) t) (-R = (^1) • ■ • 7 nv)) of an N-hodj Hamiltonian H with a particle 
density whose correct distribution is found by following the random walk 
of many test particles ("walkers") in imaginary time in 3A^-dimcnsional 
configuration space. To reduce fluctuations one works with a density distri- 
bution f{R, t) = ^{R, t)^g{R), where is a given guiding function used 
for importance sampling. The density distribution / obeys a drift-diffusion 
equation in imaginary time. Because the importance-sampled Green's func- 
tion is an exponential operator, one can expand it in terms of a Euclidean 
path integral. For sufficiently small time steps one can write down accu- 
rate approximations to the Green's function, and sample it with diffusion 
Monte-Carlo.3-6 

2.1. Choice of the guiding functions 

The choice of the guiding function is crucial for the success of the DQMC 
procedure. We take the guiding function 'J'q^ as a Slater determinant of 
single-particle orbitals each of which is a product of a Landau state in the 
lowest level with a given magnetic quantum number and an unkown longi- 
tudinal wave function ("adiabatic approximation"^). The different longitu- 
dinal wave functions are obtained selfconsistently by an iterative solution 
of the Hartree-Fock equations using B-splines on finite elements. 

2.2. Jastrow factor 

To incorporate correlation effects it is common to miiltiply the guiding 
function by a Jastrow factor, = = e"'^^-'^)*^^. We adopt the 

form 

N N 

U = -l/4Y,rij/{^ + \/Prij)+zY,n/{l + y/Pri) , (1) 

i<j i=l 

where (3 is the magnetic field strength in atomic units (/3 = B/Bq, Bq = 
4.701 X 10^ T). This leads to modifications of the adiabatic approximation 
guiding functions only at small distances, of the order of the Larmor radius. 
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3. Results cind Discussion 

As a representative example, Fig. 1 shows for the ground state of neutral 
iron (Z = 26) at B = 5x 10* T the typical flow of a diffusion quantum Monte 
Carlo simulation. Ions can be treated without additional complication in 
the same way.^ The figure depicts the energy offset E^, the block energy 
Eb and the averaged block energy (Eb) as a function of the number of 
blocks performed. 

The complete simulation goes through three stages. During the first 

100 blocks, a variational quantum Monte Carlo calculation (VQMC) is per- 
formed. Since the adiabatic approximation guiding wave function is aug- 
mented by the Jastrow factor, the VQMC calculation already lowers the 
energy in comparison with the initial adiabatic approximation result. This 
stage is followed, in the next 300 blocks, by a fixed-phase diffusion quantum 
Monte Carlo (FPDQMC) simulation. It is seen that the onset of the sim- 
ulation leads to a considerable drop in the energy. Finally, in the last 300 
blocks a releascd-phase diffusion quantum Monte Carlo (RPDQMC) simu- 
lation is carried out, which still slightly lowers the averaged block energy, by 
roughly 0.1 per cent. The dashed vertical lines in Fig. 1 indicate the blocks 
where dynamical equilibrium of the walkers is reached. The relatively small 
difference between the fixed-phase and the released-phase results indicates 
that the phase of the adiabatic approximation wave function already well 
reproduces the phase of the ground state wave function. The small fluctua- 
tions of the individual block energies Eb evident in Fig. 1 are characteristic 
of diffusion quantum Monte Carlo simulations. It is also seen, however, that 
the averaged block energies (-Eb) quickly converge to constant values in all 
three stages of the simulation. Our final RPDQMC result for the energy is 
Eq = —109.079 keV and lies well below the density fimctional (DF) value. 
The standard deviation of the block energies at the end of the simulation 
in this case is u = ±0.186 keV. 

Table 1 lists the results for all elements from helium to iron at the 
magnetic field strength B ~ 10* T. The table contains in the first three 
columns the results of the three stages of the simulation and in the fourth 
column the energy values in adiabatic approximation calculated with our 
own Hartree-Fock finite-element method HFFEM. Literature values ob- 
tained by Ivanov and Schmelcher" (2DHF), by Mori and Hailey* (MCPH^, 
multi-configurational perturbative hybrid Hartree-Hartree-Fock) and the 
results of density functional calculations^'^*^ (DF) are given in the remain- 
ing columns. The numbers in brackets attached to the HFFEM, 2DIIF, 
MCPH'^ and DF results designate the number of electrons occupying an 



Table 1. Energy values in keV for the ground states from helium to iron at B = 10* T. Parameters of the 
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Note: *Ref. 9. ''Ref. 10. 
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Fig. 1. Behavior of the block energy Eb (ragged curve) and the averaged block energy 
(E-b) (smooth curve) in the DQMC simulation for the ground state energy of neutral iron 
(Z = 26) at B = 5 X 10^ T as a function of the number of blocks. In each block, 200 time 
steps At = 5 X 10~® a.u. were performed. (HFFEM (top horizontal line); energy value in 
adiabatic approximation; DF (second horizontal line from top): density functional result 
of Ref 9; MCPH^ (third horizontal line from top): result of Ref. 8.) 



excited hydrogen-like single-particle longitudinal state. It can be seen that 
already the fixed-phase results lie slightly below the values that were ob- 
tained using the 2DHF method. The comparison with the results of the 
MCPH'^ method shows that our RPDQMC energy values generally lie be- 
low those results, but there arc also exceptions where our results lie above 
the MCPH'^ energies. This may be due to the fact that the hybrid method 
is not self-consistent, since it evaluates the exchange energy in first-order 
perturbation theory in a basis of Hartree states and it does not include the 
back-reaction of the excited Landau states whose admixtures are taken into 
account perturbatively on the effective interaction potentials. Therefore the 
method need not necessarily produce an upper bound on the energy. 

The comparison with the results of the DF calculations shows that these 
yield lower ground state energies at small nuclear charge numbers than our 
RPDQMC results, while for large Z the reverse is the case. The DF results 
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listed in Tabic 1 differ in the choice of the exchange functional. Given 
this restriction, it cannot be ensured that the DF calculations in all cases 
produce an upper bound on the ground state energy in magnetic fields as 
do the ab-initio methods used in this work or in a Ref. 11. 

4. Conclusions 

We have extendend the released-phase diffusion Monte Carlo method to the 
calculation of the ground state energies of atoms and ions from helium to 
iron neutron star magnetic field strengths by using adiabatic approximation 

wave functions as guiding wave functions.^ However, for matching observed 
thermal spectra from isolated neutron stars, wavelength information, and 
thus energies of excited states, are requisite. Jones et al.^ have shown a way 
how to calculate excited states of small atoms in strong magnetic fields using 
the correlation function Monte Carlo method.'^ The challenge remains to 
transfer their method to the DQMC simulations presented in this paper, 
and to calculate excited states of large atoms in intense fields. 
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